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Abstract 

X-ray emission from the supernova remnant N23 in the Large Magellanic Cloud (LMC) is studied 
using the X-ray Imaging Spectrometer (XIS) onboard Suzaku. Thanks to superior energy resolution of 
the XIS in the soft X-ray band, we resolved H-like and He-like Oxygen Ka emission lines from N23 with 
unprccedentedly high quality, and as a result, identified a new optically thin thermal emission component 
with a temperature —0.2 keV, as well as that with a temperature of ~0.5-0.7 keV previously known. 
This alters the estimate of the ionization timescale n e t from — 10 10-11 cm -3 s to >10 12 cm -3 s. Under 
the assumption that N23 is still in the Sedov phase, its age evaluated from the newly discovered low 
temperature component is ~8000 yr, although it is possible that N23 has already moved into the radiative 
phase. The abundances of the heavy elements are found to be roughly consistent with those of the LMC 
average, which indicates that the origin of the X-ray emission of N23 is swept-up ambient material, as 
expected from its ionization timescale. 
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1. Introduction 

Chemical evolution of the universe has been one of the 
major issues in the modern astronomy Heavy elements, 
or metals, have been generated and accumulated in the 
universe since its birth mainly through supcrnovae. An 
important clue to understand the chemical evolution of 
the universe will no doubt be brought about by system- 
atic study of a well-defined sample of supernova remnants. 
In X-ray observations, we are able to extract information 
of plasma temperature, metal abundance, explosion en- 
ergy, ionization age, and so on from optically thin ther- 
mal plasma emission from supernova remnants (SNRs). 
These pieces of information will lead to understand galac- 
tic chemical evolution and star formation history. A sys- 
tematic study of SNRs in the Large Magellanic Cloud 
(LMC) suits well for this purpose, because of its well- 
known distance (50kpc; Feast 1999) and small interstellar 
absorption to LMC. The fact that LMC is face-on to us 
is another advantage. Because of these characteristics, 
systematic X-ray studies of the thermal SNRs in LMC 
have been carried out with some major X-ray observato- 
ries (e.g. Hughes et al. 1998). These observations are, 
however, limited in terms of energy resolution in the low 
energy band below — 1 keV, which is important to study 
the nature of the SNR plasma in detail, because the tem- 
perature of the plasma in most SNRs is less than —2 keV. 
In particular, hydrogenic and He-like Oxygen Ka lines 
that appear in 0.5-0.7 keV dominate a spectrum of the 
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plasma with a temperature of <1 keV, and hence, high 
resolution spectroscopy in this energy band is of great im- 
portance in evaluating the parameters of the low temper- 
ature plasma, such as the temperature, the density, the 
ionization parameter and the abundances of the metal. 
We therefore have decided to perform a systematic study 
of the SNRs in LMC with the Suzaku XIS (Mitsuda et 
al. 2007, Koyama et al. 2007), which has the best energy 
resolution in the Oxygen 0.5-0.7 keV energy band among 
the CCDs currently in orbit. In this paper, we study the 
SNR N23, as a beginning of our systematic study project. 
N23 is identified as a SNR with radio observations at 
5 and 14.7 GHz for the first time (Milne et al. 1980). Its 
X-ray emission extends 100" x 120" in the sky (24 x 29 pc 
at a distance of 50 kpc), showing a semi-circular morphol- 
ogy with only the south-eastern hemisphere being bright 
(Hughes et al. 2006). Based on the ASCA observation, 
Hughes et al. (1998) classified N23 as a young SNR with 
an age of ~ 3800 from the spectrum of the entire remnant. 
Hughes et al. (2006) carried out spatially resolved spec- 
troscopy of N23 with Chandra, and estimated the age of 
north-western rim to be ~4600 yr. These observations 
result in consistent estimation of the temperature and 
the ionization parameter (n e t) in the bright south-eastern 
shell, which is -0.5-0.7 keV and - 10 10 -10 n cm~ 3 s, re- 
spectively. They are, however, not be able to resolve hy- 
drogenic and He-like Kcv lines from Oxygen clearly be- 
cause of limited energy resolution below ~1 keV. In ad- 
dition, Hughes et al. (2006) have revealed that there is 
a point source around the center of N23. From its 0.5- 
10 keV flux (~ 10 -14 erg cm -2 s _1 ) and power-law spec- 
trum with a photon index of 2.2, the source is probably a 
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rotation-powered pulsar and/or a pulsar wind nebula. If 
so, the progenitor of N23 is a core-collapsed massive star 
(Hayato et al. 2006). 

This paper is organized as follows. After describing 
how the observation and data reduction are carried out in 
§ 2, we report the results of spectral analysis of N23 in 
§ 3. Owing to high spectral resolution of the XIS in the 
Oxygen Ka band, we have clearly resolved hydrogenic and 
He-like Ka emission lines from Oxygen for the first time, 
which leads to identification of a new low temperature 
component. Based on these results, we discuss the plasma 
parameters and evolutionary phase of N23 in § 4. Finally, 
the summary is given in § 5. 

2. Observation and Data Reduction 

N23 was observed with Suzaku (Mitsuda et al. 2007) 
on 2005 August 16-17. The observation log is summa- 
rized in table 1. Suzaku is equipped with four modules 

Table 1. observation log. 



08:00.0 



07:00.0 



5:06:00.0 



05:00.0 



04:00.0 




Observation ID 
Observation mode 
Start Date (UT) 
Exposure (ksec) 
Count rate* (cts s _1 ) 



100003010 

Full window, no burst, no SCI 

2005/08/16-17 

7.1 

FI: 0.6 BI: 1.2 



* Count rate of FI and BI-CCD. Count rate of N23 in 
the 0.5-2.0 band within a 3'. 3 radius circle after the back- 
ground subtraction. 

of the XIS (Koyama et al. 2007). One of them adopts a 
back- illuminated CCD (BI-CCD), which is referred to as 
XIS1 having superior quantum efficiency especially in a 
low energy band, whilst the other three modules (XIS0, 
XIS2 and XIS3) utilize front-illuminated (FI) CCDs which 
have high sensitivity in a high energy band and better en- 
ergy resolution than the BI-CCD. They are adapted at the 
focal plane of the X-Ray Telescopes (XRT; Serlemitsos et 
al. 2007). They are thin- foil-nested X-ray mirrors real- 
izing a high throughput in 0.2-12 keV with a moderate 
imaging capability (1'8— 2'3 in half-power diameter). The 
Suzaku XIS image of N23 below 2 keV is shown in Fig. 1. 
Data reduction and analysis of the present data were 
carried out using the HEADAS software package version 
6.5.1. The data processed by Suzaku pipeline process- 
ing software (ver 2.0) were analyzed. We only used the 
data taken with the 5x5 editing mode in the high and 
super-high data rates because of the editing mode of 2 x 2 
has uncertainty in the response matrices of the XIS. In 
screening the data, we removed time intervals while the 
elevation angle from the night earth is less than 5°. In 
the standard data screening procedure, the time interval 
during which the elevation angle from the day earth is 
less than 20°. This results in, however, only 5.6 ksec data 
available. We thus reduce the day-earth elevation angle by 
degrees by closely comparing resulting spectra and have 
finally found that the solar X-ray contamination does not 
raise any spectral modification even if we reduce the day 



Fig. 1. Suzaku XIS X-ray image of the N23 in the bands 
0.2-2 keV. Another nearby SNR, DEM L71, is also shown. 
The images from all the four XIS modules are combined. The 
green circles with a radius of 3.'3 are integration regions of 
the source photons from N23 The background photons were 
extracted from annular region which was removed the source 
and a center regions. The calibration sources at the corners 
are masked. 

earth elevation down to 5° (DYE_ELV > 5°) for the XIS 
data. Owing to this study, the total exposure time has 
increased to 7.1 ksec. 

3. Spectral Analysis 

3.1. Extraction of Source and Background Photons 

In extracting source photons, we adopt a circular aper- 
ture with a radius of 3.'3 centered on N23 (the green circle 
centered on N23 in Fig. 1). For the background, we have 
collected photons from the entire CCD area out of the 
source extraction region, except for outer edge regions il- 
luminated by 55 Mn isotopes adapted in the camera bodies 
of all the XIS modules, the central circular region with a 
radius of ^1?5, and the 3f3 circular region centered on the 
other SNR DEM L71. We combined the spectra from the 
XIS adopting the FI-CCDs (XIS0, XIS2, and XIS3) into a 
single spectrum, and refer to it as the FI spectrum here- 
after. The BI spectrum is the same as that solely from 
the XIS1 data. 

In the spectral analysis, we use xspec version 11.3.2aj. 
The XIS response matrix (RMF) and auxiliary response 
file (ARF) are calculated using xisrmfgen (version 2007- 
05-14) and xisarfgen (version 2008-04-05), respectively 
(Ishisaki et al. 2007). The RMFs and ARFs from the FI- 
CCDs are combined with the ftools ADDRMF and mar- 

FRMF. 
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3.2. Spectral Models 

From the ASCA observations of LMC SNRs, it is found 
that the spectra of N23 is well represented by an optically 
thin thermal emission model in ionization non-equilibrium 
(Hughes et al. 1998). In characterizing the spectra of 
SNR, we thus adopt the non-equilibrium ionization model, 
or vnei model in the XSPEC model library. The vnei 
model provides the plasma temperature kT e , the ioniza- 
tion timescale n c t where t is the elapsed time since the 
shock occurs, the metal abundances, and the normaliza- 
tion parameter ^ D i J n e nndV where D is the distance 
to the source, n a and uh are the electron and hydrogen 
density respectively, and dV is the volume element of the 
SNR. 

The emission spectra are attenuated by photoelectric 
absorption due to metals contained in the interstellar mat- 
ter within our galaxy and LMC. Since the metal abun- 
dances of these two absorption components are signifi- 
cantly different, we consider them separately in the spec- 
tral evaluation as follows. We assume solar composition 
(Anders & Grevcssc 1989) for the absorbing matter in the 
Galaxy, and denote its hydrogen column density as N§ . 
The value of N§ toward N23 is obtained by the Galactic 
Hi survey, which is 5.8xl0 20 cm~ 2 (Dickey & Lockman 
1990). We utilize the PHABS model in xspec to represent 
the galactic absorption, and freeze N§ at these values in 
the spectral fitting described below. The average metal 
abundances of LMC, on the other hand, are measured 
by Russell & Dopita (1992), which is ^0.3 solar on aver- 
age over major metals. We reflect these abundances on 
the vphabs model, and set the column density associated 
with LMC iVg free to vary in the following spectral fit 
process. 

In evaluating the spectrum of N23 through spectral fit- 
ting, we begin with a single component vnei model at- 
tenuated by photoelectric absorption. Since the carbon 
abundance cannot be constrained because of the large in- 
terstellar absorption at carbon Ka line energies, we fixed 
it to the LMC value (Russell & Dopita 1992). In the early 
phase of the Suzaku mission, the energy gain has a large 
uncertainty (Koyama et al. 2007). The gain uncertainty 
is ~ ±10 eV, according to the previous studies on SNRs 
(Bamba et al. 2008; Yamaguchi et al. 2008). We therefore 
have allowed the energy offset to be floated. As a result, 
the energy offset is converged to ~ —8 eV and ~ — 1 cV, 
for the FI and BI spectra, respectively. These gain off- 
set values are within the gain uncertainty. The result 
is summarized in the second column of table 2 labelled 
TfcTini'. The best fit results in kT c ~ 0.548 keV and 
log(n c < [cm- 3 s]) ~ 10.51 with a X 2 (d.o.f.) of 144 (99). 
These parameters are roughly consistent with those from 
the Chandra observation. Note that, if we do not float the 
energy gain, the best-fit parameters are kT e ~ 0.552 keV 
and log(n e < [cm" 3 s]) ~ 10.54 with a X 2 (d.o.f.) of 219 
(101). Hence, the energy offset adjustment significantly 
improves the fit, while the differences in the tempera- 
ture and the ionization timescale (in logarithm) are within 
~1%. 




Energy (keV) 

Fig. 2. The data and the best-fit '2kTlnt' model of N23. in 
the 0.45-2.0 keV band. The black and red crosses show the 
data point from the FI and BI CCDs, respectively. The solid 
histograms with the same color are the best-fit model. The 
blue and green broken histograms show the components with 
fcT e =0.57 and 0.22 keV, respectively. The lower panel shows 
the residuals from the best-fit model. 

As noticed from this table, however, the single com- 
ponent vnei model no longer provides an acceptable fit 
(xl — 1.5). We thus have appended another vnei model. 
To reduce the number of model parameters as much as 
possible, we constrained either kT c or n c t common be- 
tween the two vnei components, and fit them to the 
data separately. The resultant best-fit parameters are also 
listed in table 2, labelled i lkT2ntf and l 2kTlnt\ respec- 
tively. Although the lkT2nt model improves the fit sig- 
nificantly, the reduced \ 2 value indicates that it is still 
not acceptable. The improvement of the 2kTlnt model, 
on the other hand, is more remarkable, and it provides 
an acceptable fit at the 90% confidence level. The best-fit 
model is shown in Fig. 2. This result indicates that there 
exists another optically thin thermal plasma component 
with a temperature of ~0.2 keV, in addition to the 0.5- 
0.6 keV component so far known. The discovery of the 
new low temperature component is brought about by the 
clear resolution of the Oxygen Ka lines by the Suzaku XIS. 
Consequently, the best-fit value of the ionization timescale 
now becomes log(n i [cm _3 s]) > 12. This is significantly 
larger than that from the ASCA and Chandra observa- 
tions (10 10 -10 u cm- 3 s). 

In the Chandra observation, more than ~ 30 % of the 
total emission comes from the eastern rim (Hughes et al. 
2006). Since our Suzaku observation covers the remnant 
entirely, the best-fit parameters we have obtained are a 
kind of mean weighted mainly on this extended bright 
rim region (see Fig. 3). 

We finally remove the constraint of the ionization 
timescale between the two vnei components and set them 
free to vary independently The results are summarized 
in the '2fcT2n£' column of table 2. This change, however, 
is not effective in improving the fit. We therefore adopt 
the 2kTlnt model for N23 for discussion in later sections. 
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Table 2. Best-fit parameters for N23* 



Parameters 


IkTlnt 


lkT2nt 


2kTlnt 


2kT2nt 


iV£(xl0 2U cm- 2 ) 


5.8(fix) 


5.8 (fix) 


5.8(ffx) 


5. 8 (fix) 


iV£(xl0 21 cm- 2 ) 


1.6(1.3-1.9) 


<0.16 


1.5(0.6-1.9) 


0.9(0.6-1.9) 


kT cl (keV) 


0.548(0.544-0.567) 


0.572(0.557-0.580) 


0.569(0.556-0.593) 


0.581(0.551-0.597) 


kT c2 (keV) 


- 


- 


0.218(0.214-0.230) 


0.225(0.214-0.230) 


log(n e i)i (cm" 3 s) 


10.51(10.48-10.55) 


10.50(10.45-10.52) 


13.68(11.92-13.70) 


12.04(11.87-13.70) 


log(n e i) 2 (cm~ 3 s) 


- 


13.18(12.77-13.70) 


- 


13.63(13.27-13.70) 


Norm^ (xl0~ 2 ) 


1.87(1.74-1.90) 


0.34(0.31-0.68) 


0.83(0.64-1.00) 


0.78(0.66-0.98) 


Norm 2 ' { (xl0~ 2 ) 


- 


0.73(0.67-0.81) 


2.54(1.93-2.96) 


1.93(1.47-3.51) 


C 


0.30(fix) 


0.30(fix) 


0.30(fix) 


0.30(fix) 


N 


<0.02 


0.09(0.04-0.25) 


0.54(0.12-1.05) 


0.37(0.23-0.94) 





0.093(0.089-0.098) 


0.27(0.26-0.28) 


0.33(0.25-0.50) 


0.34(0.28-0.49) 


Ne 


0.15(0.14-0.16) 


0.43(0.39-0.46) 


0.49(0.35-0.70) 


0.49(0.40-0.66) 


Mg 


0.14(0.12-0.17) 


0.31(0.27-0.35) 


0.39(0.33-0.60) 


0.39(0.32-0.47) 


Si 


0.28(0.19-0.36) 


0.34(0.24-0.43) 


0.41(0.30-0.52) 


0.42(0.32-0.55) 


Fe 


0.111(0.105-0.117) 


0.18(0.17-0.22) 


0.24(0.19-0.33) 


0.24(0.19-0.33) 


X 2 (d.oi.) 


1.45(99) 


1.11(97) 


1.00(97) 


1.00(96) 


offset FI (eV) 


-8.7(-11.0— 7.8) 


-8.1(-9.1— 5.9) 


-8.1(-9.0 — 7.1) 


-8.0(-9.0 — 7.1) 


offset BI (eV) 


-1.0(-1.6- 0.0) 


-1.0(-2. 1-0.0) 


-1.0(-2.8-0.0) 


-0.3(-1.3- 0.0) 



* The errors in the parentheses represent the 90% confidence intervals 

* In units of 10 — 14 /47r_D 2 / n e nndV cm~ 5 , where V and D are the Volume and distance to the plasma, respectively. 



4. Discussion 

4-1. Densities, total energy, and swept up mass 

In this section, we aim to derive electron number den- 
sities of the SNR plasma (n c ) and ambient ISM (no), to- 
tal thermal energy of the plasma {Et), swept up mass 
(M SWGp t). With the aid of the normalization obtained from 
the spectral fitting, the volume of the plasma (V) and the 
volume filling factor (/), we calculated n G , taking into ac- 
count n e =1.2riH for fully ionized solar abundance plasma. 
Assuming the strong shock, prc-shock density (e.g. ISM 
density; no) is obtained to be no = (nn/4). The swept-up 
mass is given by Af swep t = nRm p fV, where m p is the pro- 
ton mass. The total energy is given by Et = 3n c kT fV , 
where k and T is the Boltzmann constant and the tem- 
perature of the plasma, respectively. Note that we assume 
that electrons and ions temperature are in thermal equi- 
librium. The age of the plasma i ion is obtained from the 
ionization parameter (n c t) which is evaluated from the 
spectral fitting. 

In order to estimate the densities, we need to know 
the volume of the plasma from an X-ray image. Since 
the Suzaku XIS has only moderate imaging capability, 
we carried out Chandra imaging analysis first. We ana- 
lyzed the Chandra ACIS data of N23 (seq# 500212) ob- 
served on 2002 December 29-30. Images are extracted 
from a standard cleaned event file using the DMCOPY com- 
mand. Fig. 3 shows the X-ray images in the (a) 0.2- 
0.8 keV and (b) 1.5-2.0 keV energy bands. The low and 
high energy bands were chosen to show the contribution 
from the ~0.2 keV and the ^0.6 keV components, respec- 
tively. In both images, the structure of N23 shows a bright 
hemisphere with ~35" (8.5 pc) radius. We therefore as- 
sumed that the N23 plasma distributes in the shell of the 
hemisphere with a radius and a thickness of 35" and 2.7" 
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Fig. 3. Smoothed Chandra image of N23 (a) in the 
0.2-0.8 keV and (b) in the 1.5-2.0 keV band. The contours 
overlaid with labels of 1.0, 5.0, 15.0 (left) and 0.1, 0.5, 1.0 
counts pixel - 1 
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Table 3. Input and derived properties for N23. 



Table 4. The characteristic ages for N23. 



Property 



Value 



low-T component high-T component 



low-T component high-T component 
kT c (kcV) 0.22 0.57 

R (pc) 8.5 8.5 



V (10 57 cm 3 ) 


9.3 


9.3 


n c (cm -3 ) 


oi o+2.5 f -l/2 


18.0i^/oTi 1/a 


n (cm -3 ) 


7 q+0.6 f - 1/2 

'• y -i.o/o.i 


4 4 +o.3 f -i/2 

*- 4 -0. 6/0.1 


E t (10 50 erg) 


0.3±0.1fo{ 2 


0.4±0.1/ 1/ 1 2 


M SWGpt (M ) 


91 +2 f 1/2 

zl -3-/0.1 


i 9+1 f 1/2 
iz -2./0.1 


iion (yr) 


> lOOO/o 7 ! 2 


> ISOO/o 1 . 7 ! 2 



(0.7 pc), the latter of which is 1/12 of the shell radius 
expected for the ISM compressed in the case of the strong 
and adiabatic shock, which assumption is to be consid- 
ered later in § 4.2. As a result, the volume of the plasma 
becomes F = 9.3x 10 57 cm 3 . 

The other parameters are summarized in table 3. 
As shown in Fig. 3, the emission from N23 is highly 
anisotropic. We therefore normalized the volume fill- 
ing factor / by 0.1 (= /o.i)- The plasma density esti- 
mated from the low and high temperature components are 
~ 31 cm -3 and ~ 18 cm -3 , respectively, as shown in table 
3. This density estimation is roughly consistent with that 
of the southeastern rim region with the Chandra observa- 
tion (~ 10-23 cm -3 ; Hughes et al. 2006). With this nor- 
malization, the ISM density estimated from the low and 
high temperature components arc ~ 8 cm -3 and ~ 4 cm -3 , 
respectively. One may imagine that this high density 
is realized through collision of the ejecta with molecu- 
lar clouds. According to a radio observation (Banas et 
al. 1997), however, no CO emission was discovered from 
a region including N23. We thus consider that N23 ex- 
plodes into a relatively dense part of the ISM (~ 10 cm -3 ). 
It is suggested this high density is related geometrically 
with the open cluster HS 114 (Hodge & Sexton 1966) near 
the remnant (Hughes et al. 2006). Finally, total ejecta 
mass is given by Mi ow + Mhigh) where M\ ow and Mhi g h 
are the mass of low and high temperature components, 
respectively. Total ejecta mass therefore is calculated to 
be ~30 Mq. Hence the plasma chiefly comprises of the 
swept-up ISM. 

4.2. SNR phase 

In this section, we consider evolutionary phase of N23 
by comparing various timescales evaluated from observed 
parameters. We compare the Sedov time scale (isedov)i 
the plasma age (ii on ) estimated from the ionization pa- 
rameter, and the cooling time scale (i C ooi) of N23. tsedov 
is based on the concept of Sedov simple blast wave model 
(Sedov 1959), in which a supernova with explosion energy 
(Eq) expands into homogeneous ISM, which is given by 



t-Scdo 



4.3x 10 J 



R 



kT P 



i pcy V 1 keV 



-0.5 



yr. 



(1) 



^Scdov (10 4 yr) 
icooi (10 4 yr) 
*ion (10 4 yr) 



0. 
0.9/ °| 



f0.26 



0.5 
1.2/0.f 6 

> 0.2,/q.! 2 



where R is the blast wave shock radius and T c is the shock 
temperature. i coo i is defined as the time scale of the tem- 
perature distribution having started to deviate from that 
of the Sedov model, which is written as 



i CO oi = 2.7 x 10 4 



En 



10 51 erg 



"o 



1 cm 3 



-0.52 



yr (2) 



(Falle 1981). If £ coc j ~ iscdov holds, the SNR is considered 
to locate at a late stage of the Sedov phase or at an initial 
radiative phase. If, on the other hand, isodov < *cooi holds, 
the SNR is regarded as being still in the Sedov phase. We 
assume that the explosion energy is 10 51 erg in eq. (2). 

We discovered the new soft emission component with a 
temperature of ~0.22 keV, in addition to the 0.5-0.7 keV 
component so far known. Accordingly, we need to update 
the remnant age estimation. The Sedov time scale and the 
cooling time scale calculated with eq. (1) and (2) are sum- 
marized in table 4, together with the ionization age ti on 
obtained based on the spectral fit. Since the spectral fit 
with the '2fcTlni' model provides only the lower limit for 
the ionization timescale n e t, we only have a lower limit for 
iion- We thus have decided to estimate the remnant age 
with eq. (1) and (2). Using eq. (1) with R= 8.5 pc, isodov 
estimated from the low and high temperature components 
becomes ^8000 and ~5000 yr, respectively, which are con- 
sistent with the lower limit of t[ on . This is a direct con- 
sequence of eq. (1), implying isodov being proportional to 

— 1/2 

T c . The high temperature component is younger than 
the low temperature component. In the Sedov scheme, 
this can be interpreted as the high temperature compo- 
nent being heated recently. On the other hands, t coo \ of 
the low and high temperature components arc 9300/of 6 
and 12000/of 6 , respectively, using eq. (2) with n given 
in table 3 (§4.1). Since isedov ^S icooi for both components, 
N23 is likely to be at a late stage of the Sedov phase, and 
if so, its age is ^8000 yrs old from the newly discovered 
low temperature component. 

Note, however, that this discussion depends on the es- 
timation of the volume filling factor /. Accounting for 
the highly anisotropic image (Fig. 3), we have normalized 
/ by 0.1. However, if the shock wave of N23 propagates 
in a higher density ISM, like the edge of the open cluster 
HS 114 (Hodge & Sexton 1966), the total volume should 
be smaller for the given emission measure, resulting in a 
smaller /. If, for instance, / = 0.01, t coo i of the low and 
high temperature components are ^5000 and ~7000 yr, 
respectively. Since i coo i < ^Scdov, the low temperature 
component should be regarded as in the radiative phase. 
In addition, we would also like to remark that, in § 4.1, 
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Fig. 4. Elemental abundance as derived from fits to the XIS 
data. The derived N23 abundance value is the square. The 
ratios are all relative to solar (Andres & Grevesse 1989). The 
silicon abundance in Russell & Dopita's work is highly un- 
certain as indicated by our use of a large error bar for this 
species. Dashed line represents the abundance in these SNRs 
are same to the LMC average. 

we have calculated the plasma volume, assuming that the 
shell thickness is 1/12 of the shell radius, which is expected 
for a spherical strong shock. This holds only if N23 is still 
in the Sedov (adiabatic) phase. If N23 has been in the ra- 
diative phase, then the shell thickness should be smaller 
due to radiative cooling, and hence, the plasma density 
would be larger than the current estimation. As a result, 
the cooling time scale estimated from eq. (2) should be 
smaller, and i coo i ^isedov Consequently, we cannot deny 
possibility that N23 has been in the radiative phase. 

In summary, the age of N23 is estimated to be ^8000 
yr from the newly discovered low temperature component 
as long as it still stays in Sedov phase. This condition 
is, however, uncertain, because the discussion depends on 
the volume filling factor of the plasma. If, for instance, 
the shock propagate in the higher density ISM, N23 may 
have already entered into the radiative cooling phase. In 
this case, the remnant age estimated from the cooling 
timescale is smaller. 

Finally, we can make an independent age estimation 
from the central source detected by Chandra (Hayato et 
al. 2006). Its luminosity is 1 x 10 34 erg in the 0.5-10.0 kcV 
(Hayato ct al. 2006). If the central source is a pulsar (al- 
though the pulsation has not been detected), the char- 
acteristic age is ~ 10 3 -10 45 years (Possenti et al. 2002), 
which is consistent with our age estimation update. 

4- 3. Abundance 

The main component of X-ray emission originates from 
the swept-up ISM in N23. Hence the elemental abun- 
dances of N23 reflect those of its environment. Figure 4 
shows the abundances obtained from the spectral analy- 
sis (table 2). The abundances of N, O, Ne and Fe agree 
roughly with the LMC average (Russell & Dopita 1992). 
The Si and Mg abundances are smaller than those ob- 



tained (Russell & Dopita 1992), although their errors are 
relatively large. The Mg and Si abundances obtained by 
Suzaku are consistent with those from ASCA (Hughes et 
al. 1998), respectively, which seem slightly less than the 
LMC average. 

5. SUMMARY 

We observed N23 in the Large Magellanic Cloud with 
Suzaku, and discovered the new soft emission component 
with a temperature of ~ 0.22 kcV, in addition to the ~ 0.5- 
0.7 keV component so far known. This alters the esti- 
mate of the ionization parameter (n Q t) significantly from 
- 10 10 - 10 11 cm" 3 s to > 10 12 cnT 3 s. With the aid of 
Chandra imaging capability, we are able to calculate the 
density of the plasma, and have confirmed that the super- 
nova explosion of N23 occurred in a high density region 
with an ISM density of ~4-8 cm -3 . The relatively high 
ambient density may be related to the open cluster HS 
114 (Hodge & Sexton 1966) near the remnant (Hughes et 
al. 2006). The parameters of the plasma from our analy- 
sis indicate that N23 is either at a late stage of the Sedov 
phase or in the radiative cooling phase. Assuming it is 
still in the Sedov phase, we have estimated the age of 
the remnant to be ^8000 yr from the newly discovered 
soft component, which is twice as old as the estimation 
by Hughes et al. (1998). Given the uncertainty of the vol- 
ume filling factor, however, we cannot deny the possibility 
that N23 has already entered into the radiative phase. 

The abundances of N23 is roughly consistent with the 
LMC average (Russell & Dopita 1992). The Si and Mg 
abundances obtained by spectral fitting are smaller than 
those obtained by Russell & Dopita (1992), although their 
errors on Si and Mg abundances are relatively large. The 
Mg and Si abundances obtained by Suzaku are consistent 
with those from ASCA (Hughes et al. 1998). 

We thank Dr. Y. Maeda for his useful discussion and 
comments. We would like to express our gratitude to all 
members of the Suzaku team for their contributions to the 
instrument preparation, spacecraft operation, software de- 
velopment, and in-orbit instrumental calibration. We also 
thank an anonymous referee for useful comments. 
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